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For any function F(x) having a Taylor expansion we solve the boson normal ordering problem for 
F [(a^) r a s ] , with r, s positive integers, [a, a^] = 1, i.e. we provide exact and explicit expressions for 
its normal form J\f {F [(a t ) r a s ] } = F [(a t ) r a s ] , where in N (F) all a's are to the right. The solution 
involves integer sequences of numbers which, for r, s > 1, are generalizations of the conventional 
Bell and Stirling numbers whose values they assume for r = s = 1. A complete theory of such 
generalized combinatorial numbers is given including closed-form expressions (extended Dobinski 
O ■ - type formulas), recursion relations and generating functions. These last are special expectation 

Consider a function F(x) having a Taylor expansion around x — 0, i.e. F(x) — X^fcLo F fc /°^ x k . In this note we 
will collect the formulas concerning our solution of the normal ordering problem for F [(a^) r a s ] , where a, a) are the 
boson annihilation and creation operators, [a, at] = l ; and r and s are positive integers. The normally ordered form 
' of the operator F [(aJ) r a s l is denoted by Af {F [(a^) r a s ] }, where in Af(F) all the a's are to the right. It satisfies the 
CN| ■ operator identity: 

■ 

O; M{F[(alya s ]}=F[(ara s }. (1) 
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Furthermore, an auxiliary symbol : 0(a, a') : is used, which means expand O in powers of a and a? and order normally 
without taking into account the commutation rule [1] , [2] . 

For the moment we restrict ourselves to the case r > s, the other alternative being treated later. We do not give the 
proofs of the formulas here; they will be given elsewhere [3] . The case r = s = 1 is known [1] , [2] and some features 
of the r > 1, s = 1 case have been published [4]. When needed we shall refer to [1], [2] and [4] where particular cases 
of our more general formulas are discussed. 

We define the set of positive integers S r ^ s (n, k) entering the expansion: 

cr. 



[(o t ) r o a ]" = AA{[(a t ) r a s ]"} = ( a t)™(^) 



(2) 



for n = 1,2, ... . Once the S rtS (n,k) are known, the normal ordering of [(a^) r a s ]™ is achieved. The same applies to 
any operator F [(a^) r a s ~\ if F{x) has a Taylor expansion around x = 0. The row sums of the triangle S r . s (n,k) are 
given by: 

ns 

B r , s (n) = ^ S r,s(n, k), (3) 

k—s 

which for any t extend to the polynomials of order ns defined by: 

ns 

B r , s {n,t) = J2 S rA^k)t k . (4) 
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All the subsequent formulas are consequences of the following relation, linking the polynomial of Eq.(4) to a certain 
infinite scries: 



oo n 

^Efcill^+y-^-^-^+O'- 1 )^-*)- 1 ) • 

...■(k+(j-l)(r-s)-s + l)]t k = 

ns 

= J2S r , s (n,k)t k , n = l,2... 



A -s ' j = 



(5) 



which for r, s > 1 and t = 1 is an analogue of the celebrated Dobinski relation, which expresses the combinatorial Bell 
numbers #1,1(77) as a sum of an infinite series [5], [6], [7]: 



oo 



A:=0 



k n 



(6) 



For a review of various characteristics of Bell numbers -81,1(71) see [8]. It seems that Bell and Stirling numbers are 
now beginning to appear in textbooks of mathematical physics [9] . 

By setting t — 1 in Eq.(5) we obtain equivalent forms of expressions for the generalized Bell numbers B rs {n),n = 
1,2...: 

r > s : B rtS (n)= B r ^ s (n, 1) = 
1 ~ 1 



;Ey l[[(k+(j-l)(r-s)).(k+U-l)(r- 3 )-l) 

...■(k + (j-l)(r-s)-s + l)] 



k=s j = l 

1^ 1 n ^ (k+jr-(j-l)s)\ 



(r - s) 



r- 

s(n-l) 



y- 

^ k\ 



k=0 



A£(n+J±f) 

y r(i + m) 



1 00 1 

= s: B rtr (n)= B r>r (n,l) = - E p 



e — jfc! 

A;=0 



(fc + r)! 



fc! 



(7) 
(8) 

(9) 
(10) 

(11) 



For all r, s we set B r ^ s (0) = 1 by convention. In Eq.(8) we have employed the notation m- := m ■ (m — 1) • (m — 2) • 
. . . • (m — s + 1) for the falling factorial [10], and in Eq.(10) T(y) is the Euler gamma function. 

The numbers S r ^ s (n, k) are non-zero for s < k < ns, and satisfy by convention S r , s (n, 0) = £„ ; o- We shall refer to 
them as generalized Stirling numbers of the second kind. Their exact expressions in the form of a finite sum are: 



S r , s (n,k) = K ) T7 [(p+ (j - i) (r _ s) ) . (p+ (j - i)( r - S ) - i) . 



...■(p+(j -!)(»■ -«)-« + !)] 
JJ(p+0--l)(r-s))* 



(-1)' 
jfc! 



p=0 



(12) 
(13) 

(14) 



x=l 



from which for r = s = 1 one obtains the standard form of the classical Stirling numbers of the second kind [7]: 
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^,i(n,fc) = ^E(- 1 ) p (p)^ ( 15 ) 
p=i ^ ' 



and for r = 2, s = 1, 



52 > l(n ' fe) = ^(fc-l)' 



(16) 



which are the so-called unsigned Lah numbers [4], [7]. 

That the numbers S ryS (n, k) are natural extensions of Si t i(n, k) can be neatly seen by observing their action on the 
space of polynomials generated by falling factorials. Let x be any real number. Then, 

n 

U(x + (j - l)(r - *)) • (x + (j - l)(r - s) - 1) • . . . • (x + (j - l)(r -s)-s + l) = 



= SrA n , k ) x ( x - 1) ■ • • (a; - fc + 1), (17) 

k—s 

which, when rewritten with m- symbols become: 

n ns 

]J(x + (j - l)(r -*))* = £ S r , a (n, k)xK (18) 

J = l fe = s 

This last equation when specified to r = s gives a particularly transparent interpretation of S r ^ r (n, k) as connection 
coefficients: 

nr 

= s r , r (n, fc)a£. (19) 

By choosing r = s = 1 it boils down to the defining equations for Si i(n, fe) in terms of x-, used by J. Stirling himself 
[10]: 



x n = ^S hl {n,k)xK (20) 
fe=i 

The relations for S r ^ s (n 7 k) become particularly appealing for r — s: Eq.(13) can be reformulated as: 

SrAn, k) = J2(-l) p (*)|p(p-l)...(p-r + 1)]" (21) 



= ^X> 1 ) P (jiW, (22) 

p— r 

which differs from Eq.(15) in that the sum in Eq.(22) starts from p = r and p n becomes [p 1 -]™. The recurence relations 
for S r , r (n, k) are the following: 

5^(1,^ = 1, S rtr {n, k) = 0, fc<r, nr < fc < (n + l)r, (23) 

5„(n + l,fc) = ^](' fc+P_r ")rP5' r!r (n,fc+p-r), r < k < nr, n > 1, (24) 

where the definition r- := 1 is made. Note that for r = s = 1 we get the known relation for conventional Stirling 
numbers i.e.: Si ; i(n + 1, k) = kS\^(n, k) + k — 1) with appropriate initial conditions [10]. 

The "non-diagonal" generalized Bell numbers B r s (n) can always be expressed as special values of generalized 
hypergeometric functions p F q . Algebraic manipulation of Eqs.(7)-(11) yields the following examples: 

r > l,s = 1: 
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B r ,i(n) is a combination of r — 1 different hypergeometric functions of type iF r _i(. . . ; x), each of them evaluated 
at the same value of argument x = (r — l) 1_r ; here are some lowest order cases: 



B 2,i(n) 
B 3 ,i(n) 

B 4 ,i(n) 



1 f 1 (n+l;2;l) = (n-l)!Li 1 2 1 (-l), 



2I> + i) 



113 1 



2' 2' 2' 4 



1^2 n+w.^T +"!i^2 n+l;-,2;- 



3 „ 1 



2' '4 



(25) 
(26) 



gn-l / 3 3/2 r( 2 )r(n+ l ) 

2e I 

3r(n+|) 



1^3 \n + -, 



112 4 1 

3' 3' 3' 3' 27 



+ 



r( 



2 2 4 5 1 



(27) 



etc. In Eq.(25) Lm\y) is the associated Laguerre polynomial. 

Similarly the series B2 r , r (n) can be written down in a compact form using the confluent hypergeometric function 
of Kummcr: 



B 2r , r (n) 



(rn)\ 
e ■ rl 



1 F 1 (rn + l,r+ 1;1). 



(28) 



A still more general family of sequences arising from Eqs.(7)-(11) has the form (p,r = 1, 2 . . .): 



pr-\-p,pr 



(n) 



n 



(p(n-l)+j)! 
(Pi)! 



r F r (pn + l,pn + l+p,. .. ,pn + 1 +p(r — 1); 1 +p, 1 + 2p. . . , 1 + rp; 1), 



(29) 



etc. 



In contrast, the "diagonal" numbers B r ^ r (n) of Eq.(ll), which also can be rewritten as: 

1 oo 

Brr(n) = _ £ _____ _ [k[k + l)...( k + r- 1)]" , n = 1, 2, . . . 



(30) 



cannot be expressed through hypergeometric functions. However, the sequences B. r _ r (n) possess a particularity that 
they can be always expressed in terms of conventional Bell numbers with r-nomial (binomial, trinomial. . . ) coefficients. 
For example: 



i? 2 , 2 (n) = E( n fc 1 ) B M(" + fc )- 



(31) 



Several low-order triangles of S r ^ s (n,k) and their associated B r _ s (n) are presented in Table 1. 

It turns out that various generating functions for B r ^ s {n) and S r ,s{n,k) can be related to special quantum states, 
called coherent states [1], defined as linear combinations of the eigenstates of the harmonic oscillator, H = a^a, H\n) = 
n\n), (n\n') = <5„ >n / and defined as: 



| z |2 ^ gll 



n=0 



(with (z\z) = 1), for z complex. The states \z) satisfy: 

a\z) = z\z). 

It has been noticed in [12] that, 

(z\{cta) n \z) 't^yfn) 



(32) 



(33) 



(34) 



and the exponential generating function (egf) of the numbers B\ t \{n), i.e. _^Lo -^i,i( n ) s& ti snes 
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{z\e^ a \s) = (z\ : e °M^-D . \ z) I'Jr 1 e e*-i = £ B M („)* 

71. 

n=0 



which is a restatement of the known fact that [1], [2], [4]: 

A/"(e Aata ) =: e a Me A -i) : 



How does it generalize to r, s > 1 ? 
Eqs.(2) and (3) give directly: 



(2|[(atyV]>>*=iMn), 



(35) 



(36) 



(37) 



valid for all r, s. 

We first treat the case r = s = 1,2.... Note in this context that a Hermitian Hamiltonian (a^) r a r is of great 
importance in quantum optics, as for r = 2, 3, . . . it provides a description of non-linear Kerr-type media [11]. Using 
Eq.(22) we can find the following egf of S r , r (n, k): 



n—\k/r~\ 



( *) fc ^(-1) P ( h \ (eMp-V-ip-r+i) _ i 



(38) 



where \y \ is the ceiling function [10], defined as the nearest integer greater or equal to y. This yields via exchange of 
order of summation: 



a Ap(p-l)...(p-r+l) 



-') 



(a T a) 



and 



(~l) fc 

fc! 



(39) 



(40) 



Another case which can be written in a closed form is r > 1, 5 = 1, for which the egf for j5 r) i(n, k) reads: 

00 n 1 

E ^^,^ = -{(1-^-1^-^-1}. 



n—\kjr~\ 

One then obtains [4], [13]: 

TV (e A ( Qt ) ra ) =: exp { [(1 - A(a t ) r " 1 (r - 1))"- 
For arbitrary r > s the following formula is still valid: 

A™ 



~ — 1 



t.}: 



N (e^^) = 1+ : ^ ^(at)^-) ]T S r , s (n, fc)( a ta)* 



(41) 



(42) 



(43) 



n=l 



\k—s 



where the explicit form of S r , s (n, k) may be used, see Eqs.(12)-(14). 
A close look at Eqs.(40) and (37) reveals that 

n=0 



(44) 



Comparing Eqs.(44),(42) and (40) we conclude that for r > 1, s = 1 and r = s the egf of respective generalized Bell 
numbers are special matrix elements of e A ( at ^ aS in coherent states. 

However, for arbitrary r > s > 1 examination of Eqs.(7)-(10) confirms that the numbers B ryS {n) increase so rapidly 
with n that one cannot define their egf's meaningfully. In particular, this is reflected by Eq.(43) which is true in its 
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operator form, but one needs to exercise care when calculating its matrix elements as the convergence of the result 
may require limitations on A. 

A well-defined and convergent procedure for such sequences is to consider what we call hypergeometric generating 
functions (hgf), i.e. the egf for the ratios B r;S /(n!)*, where t is an appropriately chosen integer. A case in point is 
the series #3,2 (^) which may be written explicitly from Eq.(29) as: 



B 3 , 2 (n) = -£ 



1 (n + k)\{n + k + 1)! 
e h / o k\(k + l)\(k + 2)\ ' 



(45) 



n=0 



#3,2 (n) 



n! 



Its hgf G 3:2 (A) is then: 

oo 

G 3 , 2 (A) = ]T 

n 

Similarly for G^2 (A) one obtains: 

1 ^ 1 

Eq.(28) implies more generally: 



\n 1 I 



(46) 



— ,- + l;l;4A). 



(47) 



n=0 



B 2r , r (n) 



(n 



W-i 



j_Y^oo 1 p f k+l fc+2 fc+3-1 1-97\A o 

2! e LtO (t+2)! 3 2 l 3 ' 3 ' 3 ' 1 -> 1 ' ZIA )' ' ~ °> 

3 iEr=0(feT3)T4F3(^l,...,^;l ) l,l;256A), r = 4...etc. 



(48) 



See [14] for other instances where this type of hgf appear. 

All the previous considerations also apply to the case r < s. In this case we define the generalized Stirling numbers 
by (note the difference with Eq.(2)): 



r < s : 



J2SrA n > k )(^) k a k 



.k—r 



a 



n(s—r) 



(49) 



By taking the Hermitian conjugate of Eq.(2), with the change r <-> s, we obtain 



r < s : 



k—r 



(1 



n(s—r) 



whence the symmetry relation: 

Sr,s{ n > k) = S s ^ r (n, k), 



s < k < ns, (r < s). 



(50) 



(51) 



The normal ordering of F[a s (a' f ) r } involves the antinormal to normal order, in the terminology of [15]. With this 
in mind we define the anti-Stirling numbers of the second kind S r ^ s (n, k) as: 



r > s : 



[aV)T = (a f r (r ~ s) 



^S s , r (n,fc)(a t ) fe a fc 



,fc=o 



It can be demonstrated [3] that they satisfy the following symmetry properties: 



&s,r(w 7 k^ — *S r , s (u, &), 

S r ^ s (n, k) = S r , s (n + 1, k + s), 



< k < ns, (r > s), 
< k < ns, (r > s). 



(52) 



(53) 
(54) 
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We have obtained recursion relations for B rjS (n) for certain values of r, s. They include the r > 1, s = 1 case: 



n / \ n — k 




(55) 



fe=0 v 7 j=0 



which for r = 2, s = 1 may be written explicitly as: 



£ 2 ,i(«+l) = E( )(n-fc + l)!B 2il (A:) 



fe=0 




(56) 



and for r = s = 1 reduces to the known relation [5], [8], which is the binomial transform: 



B M (n+l) = ]T I )B hl (k) 



k=0 




(57) 



If in Eq.(l) F(x) has a Taylor expansion of F(x) around xo ^ 0, then 



Af{F[(a^a s ]}=Af(Y t 




F^(x D ) 



(58) 



and our above results should be supplemented by the appropriate binomial expansion coefficients. 

All the B r:8 (n) described here are the n-th moments of positive functions on the positive half axis; some solutions 
of the associated Sticltjes moment problem have been recently discussed at some length [16]. In addition we have 
found that even larger classes of combinatorial sequences are solutions of the moment problem and they can be used 
for the construction of new families of coherent states [17]. 

We are considering the possible combinatorial interpretation of the above results. 
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Table 1. 

Triangles of generalized Stirling and Bell numbers, as denned by Eqs.(2) and(3). 
r=l, s=l 







Si 


i(n, k), 1 < k < n 


Bi,i( 


n = 


1 


1 




1 


n = 


2 


1 


1 


2 


n = 


3 


1 


3 1 


5 


n = 


4 


1 


7 6 1 


15 


n = 


5 


1 


15 25 10 1 


52 


n = 


6 


1 


31 90 65 15 1 


203 



n) 



r=2, s=l 

5 2 ,i(n,fc), 1 < k < n B 2 ,i{n) 



n 


= 1 


1 










1 


n 


= 2 


2 


1 








3 


n 


= 3 


6 


6 


1 






13 


n 


= 4 


24 


36 


12 


1 




73 


n 


= 5 


120 


240 


120 


20 


1 


501 


n 


= 6 


720 


1800 


1200 


300 


30 1 


4051 



r=3, s=l 



r=2, s=2 









S3,i(n, k), 1 < k < n 


£3,1 (n) 








n = 1 


1 




1 








n = 2 


3 


1 


4 








n = 3 


15 


9 1 


25 








n = 4 


105 


87 18 1 


211 








n = 5 


945 


975 285 30 1 


2236 








n = 6 


10395 


12645 4680 705 45 1 


28471 








52,2(71, fc), 2 < fc < 2n 






S 2 , 2 (n) 


n — 


1 


1 








1 


n = 


2 


2 4 1 








7 


n = 


3 


4 32 38 


12 


1 




87 


n = 


4 


8 208 625 


576 


188 24 1 




1657 


n = 


5 


16 1280 9080 


16944 


12052 3840 580 40 


1 


43833 


n = 


6 


32 7744 116656 


412800 


540080 322848 98292 16000 


1390 60 1 


1515903 



r=3, s=2 







S3,2(n, k), 


2 < k < In 














£3,2(71) 


n 


= 1 


1 


















1 


n 


= 2 


6 


6 


1 














13 


n 


= 3 


72 


168 


96 


18 


1 










355 


n 


= 4 


1440 


5760 


6120 


2520 


456 


36 


1 






16333 


n 


= 5 


43200 


259200 


424800 


285120 


92520 


15600 


1380 


60 


1 


1121881 


n 


= 6 


1814400 


15120000 


34776000 


33566400 


16304400 


4379760 


682200 


62400 


3270 90 1 


106708921 
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r=3, s=3 

53,3 (n, k), 3 < k < 3n B 3 ^(n) 



n 


= 1 


1 




















1 


n 


= 2 


6 


18 


9 


1 














34 


n 


= 3 


36 


540 


1242 


882 


243 


27 


1 








2971 


n 


= 4 


216 


13608 


94284 


186876 


149580 


56808 


11025 


1107 


54 


1 


513559 


n 


= 5 


1296 


330480 


6148872 


28245672 


49658508 


41392620 


18428400 


4691412 


706833 


63375 3285 90 1 


149670844 
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